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Abstract 



In a low-temperature study with a scanning tunneling micro- 
scope (STM), the irreducible lateral motion of a CO molecule ad- 
sorbed on a Si(OOl) surface showed a hyperlinear dependence on 
the tunneling current. This dependence implies that the adsor- 
bate displacement is caused by multiple excitations of adsorbate 
vibration modes, a situation thus far observed only at metal sur- 
faces. The local vibronic temperature at the atomic scale on the 
surface heated by ohmic inelastic scattering of tunneling electrons 
indicates that there is an activation barrier of 0.11 eV for the irre- 
versible motion of CO, in agreement with the adiabatic potential 
obtained from first-principles calculation. The highly efficient local 
heating is caused by a mid-gap state at the surface induced by the 
electric field of the STM tip. 



1 Introduction 



Techniques of atom-by-atom manipulation, including desorption, hopping, 
and chemical reaction, to modify surfaces at nanometer scale with scan- 
ning tunneling microscopes (STM) are well-established [H 121 |3]. High 
electron density of the tunneling current plays an important role in these 
surface reactions, which are typically interpreted in terms of two mech- 
anisms, (i) In the Menzel-Gomer-Redhead (MGR) model, an electronic 
excitation to a stable anti-bonding state leads to nuclear motion along 
an adiabatic potential [U [5] . This usually occurs when the applied bias 
voltage is relatively high [B], and the reaction rate R shows a linear de- 
pendence on the tunneling current / [TJ (HJ |9] . (ii) In the ladder- climbing 
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model [To], energy quanta of a local vibration are multiply excited by the 
inelastic scattering of tunneling electrons applied at a relatively low bias 
voltage, until the accumulated energy exceeds the activation energy. In 
this case, the reaction is governed by the power-law R ~ where n 
is the reaction order. Such multiply excited process have been observed 
chiefly on metal surfaces thus far [II1[I2], but for very limited cases [3|[T3]. 
Accumulation of vibrational energy at molecules on surfaces is competi- 
tive against dissipation into equilibrium. For instance, the lifetime of the 
C-0 stretching vibration on Pt(lll) is 2.2 ps [H], whereas that of CO 
on Si(OOl) is 1.87 ns [T51 [TB]. The latter longer lifetime suggests that the 
accumulated energy would overcome the reaction barrer, but the observed 
rate in vibration-induced reaction follows a linear relation on a semicon- 
ductor [I7j. Because the mechanism is still unclear, we focus in this work 
on the role of covalent electronic states localized near the molecule ad- 
sorbed on the surface. Our detailed STM experiments combined with 
first-principles calculations reveal that the accumulated energy of CO vi- 
brations on Si(OOl) during STM observation is sufficient to induce induce 
a change of the adsorption structure of the CO molecules. We will show 
that the power law breaks down in this system, and that the vibrational 
temperature is a good indicator of atomic-scale reactions proceeding on 
the semiconductor surface. 

The Si(001)-c(4x2) surface adsorbs CO molecules below ~200 K with- 
out decomposition [181 [IH] ■ There are two stable adsorption sites [201 Ell 
|22] . where CO either (i) terminates a dangling bond at the down-Si of the 
dimer (T-CO), or (ii) sits on the metastable bridge site of two Si atoms 
of a dimer (B-CO). On a clean surface, the formation of B-CO requires 
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an activation energy of 1.1 eV, which can be achieved, for example, by 
an accelerated CO molecular beam ^23j. In our STM images of the CO- 
adsorbed surface, we show that bright spots emerge after scanning, which 
indicates that B-CO is formed in response to the tunneling current. We 
discuss this process in terms of an inelastic tunneling mechanism, based 
on density functional calculations. 



2 Experimental and Computational Details 

We used on an ultrahigh vacuum (UHV) chamber with a base pressure of 
less than 5 x 10^^ Pa, equipped with an STM (Omicron VT-STM). Our 
STM tip was made by electrochemical etching of a polycrystalline tungsten 
wire and further cleaned by electron bombardment. Before CO exposure, 
a clean c(4x2) structure of Si(OOl) (n-type, 0.02 fi-cm) was confirmed by 
STM observation, and the sample was cooled to ~90 K. The clean surface 
was exposed to CO at a level of more than 100 Langmuir |25], and STM 
observation was performed at various tunneling currents. To model the 
CO adsorption, density functional calculations were performed with the 
VASP code [261 122] • The energy cutoff of the plane wave basis set was 
set at 400 eV. We modeled the surface of a (2x3) unit cell as containing 
eight layers of Si separated by more than 10 A of vacuum. CO molecules 
are adsorbed at the reconstructed top layer, while the dangling bonds in 
the unreconstructed bottom layer are terminated with hydrogen atoms. 
For the Brillouin-zone integration, a 4x4 k-point grid was used [28]. All 
atoms in the top seven layers were relaxed until the force on each atom 
fell below 0.02 eV/A. 
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3 Results 



Fig. 1 shows successive STM images of a 50x50 nm^ area after the CO 
exposure was stopped. A few bright spots and darker islands were ob- 
served, as reported previously [21] . We found that the number of bright 
spots increased upon repeated scanning even after the CO exposure was 
stopped. The seventh scan covered a wider area of 100x100 nm^, as shown 
in Fig. 1(e). Here, the number of bright spots is obviously increased in 
the marked area, which had already been scanned six times, indicating 
that formation of the bright spots was associated with the repeated STM 
scanning. The bright spots did not appear on dimers where the adjacent 
dimer showed a bright spot. Disappearance of the bright spots was rarely 
observed. 

The magnified STM image in Fig. 2(a) shows that the bright spots 
are formed at the centers of the Si dimers, and thus can be interpreted as 
B-CO. On the other hand, the T-CO structure was not clearly observed, 
suggesting that T-CO is invisible to STM at the applied bias, because the 
shape of the electron distribution is similar to that of a dangling bond 
of intact Si dimer. A schematic structure of one B-CO and two T-CO 
in a 2x3 unit cell is superimposed on images of the electronic isosurface 
mapping calculated using the Tersoff-Hamann model in Fig. 2(b) and (c). 
Fig. 2(b) was calculated for a charge neutral system (A^e), and Fig. 2(c) 
for a system with one excess electron injected (A'e + 1) due to the electric 
field effect induced by the bias between the sample and the STM tip [29] . 
B-CO appears as a hemisphere in the experimental image. The calculated 
image of A'e clearly splits into the two lobes of a p-orbital, while that of 
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A'^e + 1 is much closer to a hemisphere. Among calculated systems with B- 
CO, only the configuration in which one B-CO is sandwiched by two T-CO 
pointing in opposite directions was stable after structural optimization. In 
all other cases, B-CO moved to the terminal site during the optimization 
due to the lack of an activation barrier against transformation from B- 
CO to T-CO. This is presumably why no adjacent B-COs were found in 
STM images after repeated scans. The energy of T-CO is lower than that 
of B-CO, which is consistent with the fact that the number of B-CO is 
negligible in the experimental image of the first scan after CO exposure. 
Thus, the emerging bright spots indicate a lateral displacement of CO, i.e., 
T-CO is initially formed during the gas exposure, and then transformed 
to B-CO due to the STM scanning. 

The reaction rate [R) of displacements is plotted as a function of the 
tunneling current in Fig. 3. Because the rate is not proportional to /, the 
reaction cannot be explained in terms of the MGR model caused by a single 
electron transition. This suggests that the reaction is a multiple excitation 
process. However, as shown in the inset, we could not fit our experimental 
data with a simple power law {R oc J"). Therefore, we introduced a local 
vibrational temperature Ty to describe the degree of vibration excitation 
so that the reaction rate can be expressed by an Arrhenius-like relation 

m 

R = Aexp{-^). (1) 

Here h is Planck's constant, Ea is the activation energy for the reaction, 
kB is Boltzmann's constant, and the pre-exponential factor A is related to 
the vibrational excitation rate and the number of multiple excitations of 
vibrational quanta spaced by hoj. This equation leads to a simple power 



law in the case where the molecular orbitals have a large overlap with sub- 
strate states near the Fermi level on metallic surfaces [SUIEI]- In contrast 
to this frequently employed approximation, energy dissipation of tunnel- 
ing electrons is described as P = 7/^ [32j , where 7 is the ohmic coefficient 
of local vibrational heating due to resistance at the tunnel junction. The 
vibrational temperature is given by % = Tq + jP, using the substrate 
temperature Tq. Hence we obtain 

R = Ae~ "Bi^olii-^) ^ (2) 

which excellently fits the curve in the main panel of Fig. 3. Note that A 
and 7 contain the thermal dissipation of energy and the frequency of in- 
elastic scattering, respectively. The obtained parameters are A = 10'^~10^ 
s-i, 7 = 14.3 ± 7.3 K/A^, and Ea = 0.11 ± 0.05 eV. The activation en- 
ergy Ea is associated with the molecular motion of CO. The calculated 
adiabatic potentials for the motion of CO along the [110] direction (Fig. 
4(a)) for no additional charge {N^) and one additional electron {Nf, + 1) 
are shown in Fig. 4(b) and (c), respectively. The calculated activation 
energy for the transformation from T-CO to B-CO is reduced from 0.17 
to 0.15 eV when the tip-induced charging of the surface is considered, in 
good agreement with the experimental results. 



4 Discussion 



Due to multiple excitations of local vibrations, T-CO transforms to B-CO 
with a potential barrier of 0.154 eV, much smaller than the value previ- 
ously calculated in a small cluster [23], indicating that repulsion between 
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CO molecules on the surface (included in our model) is non-negligible. In 
fact, the barrier is small enough that three or four excitations of vibra- 
tional quanta at T-CO overcome it. The reaction rate R obtained with 
STM does not fit with J" for any n, but is well explained by vibrational 
heating proportional to the square of tunneling current /. This result also 
indicates that an increase of temperature by several tens of K is sufficient 
to induce this structural change. Moreover, high-resolution electron en- 
ergy loss spectroscopy shows no peak related to B-CO at ~80 K [HI [19], 
while at higher temperature a shoulder of B-CO [23j appears on a T-CO 
peak [ISl This experiment suggests that the formation of B-CO can 
be induced at thermal equilibrium, but the STM heating is ultimately 
localized at the adsorbed site [311 112] • As there are many vibronic modes 
at the sites, the vibrational temperature is a good indicator for reactions. 
Indeed, this idea is widely accepted for reactions by electronic excitation 
on metal surfaces [3T]. In this inequilibrium model, the temperatures of 
electrons, the local vibration of adsorbate and the substrate phonons af- 
ter electronic excitation by photons reach equilibrium separately [35j . Our 
results show that introducing the vibrational temperature is also valid for 
surface reactions induced by STM current [36]. In our semiconductor case, 
the temperatures of CO and Si substrate are different under local current 
heating with the STM, and the efficiency of this heating depends on the 
local structure at the atomic scale. On a semiconductor surface, the vibra- 
tional energy dissipates into the bulk, and the dissipation is dependent on 
the vibrational modes [37]. The lifetime of the T-CO vibration is as long 
as 2.3 ns [15], which is comparable to the average period of electron/hole 
injection in our experiment (an electron per ns, approximately). 
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The transformation from T-CO to B-CO is irreversible, which suggests 
that the temperature is not homogeneous in an area wider than the STM 
resolution, as homogeneous heating would lead to a constant ratio of T- 
CO and B-CO. An STM tip at a bias voltage of 1.6 eV introduces 1.3 
electrons in a unit cell on Si (electronic susceptibility of e ~ 12eo) [25] . 
corresponding to our calculation of the A^e + 1 state. We find a surface 
localized state in the band gap for the T-CO structure. On the other hand, 
the extra electron in B-CO delocalizes to the bulk. A surface localized 
state is known to make a large contribution to inelastic tunneling in the 
case of the yr-orbital in benzene adsorbed on Si(OOl) [SH], and thus we 
expect that heating at the T-CO site would occur more efficiently than at 
the B-CO site. 

The power law of surface chemical reactions has been successfully 
applied for metal surfaces, where excitation only within the adsorbed 
molecule is considered. Dissipation of energy at the metallic surfaces is 
very rapid, while the lifetime of molecular excitation is much longer. In 
the STM-induced process on semiconductor surfaces, however, the excita- 
tion at the adsorbate and dissipation into the substrate are competitive. 
The vibrational energy of B-CO, obtained in our calculation, is at reso- 
nance with the Si-bulk phonon continuum, whereas that of T-CO is off 
the resonance [39j. Therefore, the vibrational energy of B-CO is strongly 
coupled to the bulk phonons, and decays very rapidly. The probability for 
transformation from B-CO to T-CO is expected to be extremely low due 
to the short lifetime of the B-CO vibration. 

The mechanism of electron transfer is described in terms of hole injec- 
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tion, as shown schematically in Fig. 5. The inelastic tunneling process at 
T-CO is very efficient because of the mid-gap state locally induced by the 
electric field E from the STM tip. The potential near CO is changed by 
both the hole injection into the local state and its recombination with a 
donor state that follows. During these transfer processes, energy of the or- 
der of ~ 150 meV can be lost, which corresponds to the vibronic quantum. 
This inelastic scattering resonantly excites the local vibration involved 
in the CO movement. At B-CO, in contrast, the hole transfers into the 
substrate, where it immediately diffuses deep into the bulk without an 
inelastic process at the surface. 



5 Summary 



T-CO adsorbed on a Si(OOl) surface transforms to B-CO during STM 
experiments due to electronic excitations. The rate of this reaction does 
not follow a power law, but refiects multiple excitations of local vibration 
modes. The local temperature, an indicator of locally accumulated energy 
at the atomic scale, follows the ohmic relation (x P. In this process, the 
key feature is vibronic excitation by inelastic scattering of an electron/hole 
tunneling through the local state induced by the tip field at T-CO. The 
long lifetime of the CO vibration allows multiple vibronic excitations to 
occur. In a semiconductor system, we have thus demonstrated that ad- 
sorbate structure modification can be induced by atomic-scale heating 
effectively owing to the STM tip. 
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Figure 1: (a-c) STM images obtained in the first to third scans of a CO- 
adsorbed Si(OOl) surface after CO exposure. 125 s was required to obtain 
each image. The scanned area is 50x50 nm^. (d) An STM image obtained 
at the seventh scan after cessation of CO exposure, (e) An STM image 
of a wider area (100x100 nm^) taken after the seventh scan to show the 
surrounding region. All images were obtained with a sample bias voltage 
of -1.6 V and tunneling current of 200 pA. 
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(b) 



(c) 



Figure 2: (a) An STM image showing CO adsorbed on Si(OOl) at -1.6 
V. (b) and (c) show the simulated structures of T-CO and B-CO with 
and A'^g + 1 electrons, respectively. Filled and open circles represent CO 
molecules and Si atoms, respectively. Below: the simulated images with 
the corresponding schematic superimposed. 




Tunneling Current (nA) 

Figure 3: Frequency of emerging bright spots in STM images during the 
average residence time of the tip at a dimer. The tunneling current was 
set to 100, 300, 600, 900 pA or 1.2 nA. Inset is the same data but plotted 
in log-log scale to show that the points do not lie on a straight line (See 
text). 
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Figure 4: (a) Schematic model showing the transition from T-CO to B- 
CO on Si(OOl). The calculated potential curves on (b) the neutral slab Ne 
and (c) the charged slab A'^e + 1 have different barriers. The calculation 
was done for a (2x3) unit cell with two fixed B-TO and one mobile. The 
reaction coordinate on the horizontal axis corresponds to the position of 
CO. The calculation involves an error of a few meV because the bulk 
atoms are not allowed to relax. 
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recombination 




Figure 5: Top: schematic band diagrams for T-CO (left) and B-CO (right) 
in the case of negative sample bias voltage. The dashed lines among bands 
indicate the sample Fermi level [Ep). Bottom: models of the impact of 
electron transfer as hole injection. The hole tunnels to the mid-gap surface 
state near T-CO then goes into the bulk region, while at the B-CO site 
the hole tunnels to the substrate without being trapped near the surface 
region. 
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